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Abstract. Many real networks are embedded in space, where in some of them the 
links length decay as a power law distribution with distance. Indications that such 
systems can be characterized by the concept of dimension were found recently. Here, 
we present further support for this claim, based on extensive numerical simulations 
for model networks embedded on lattices of dimensions d e = 1 and d e = 2. We 
evaluate the dimension d from the power law scaling of (a) the mass of the network 
with the Euclidean radius r and (b) the probability of return to the origin with the 
distance r travelled by the random walker. Both approaches yield the same dimension. 
For networks with 5 < d ei d is infinity, while for 6 > 2d e , d obtains the value of 
the embedding dimension d e . In the intermediate regime of interest d e < S < 2d e , 
our numerical results suggest that d decreases continously from d — oo to d e , with 
d—d e ~ (S — de)^ 1 for S close to d e . Finally, we discuss the scaling of the mass M and 
the Euclidean distance r with the topological distance £ (minimum number of links 
between two sites in the network) . Our results suggest that in the intermediate regime 
d e < 8 < 2d e , M(£) and r{£) do not increase with £ as a power law but with a stretched 
exponential, M(l) ~ exp[A£ s '^- s '^} and r{£) ~ exp[B£ s '( 2 - s \ where 5' = 5/d e . The 
parameters A and B are related to d by d = A/B, such that M(£) ~ r(£) d . For 
S < d e , M increases exponentially with £, as known for 5 = 0, while r is constant 
and independent of £. For 6 > 2d e , we find power law scaling, M(£) ~ £ dt and 
r{£) ~ with di ■ d m in = d. For networks embedded in d e = 1, we find the 

expected result, df = d m i n = 1, while for networks embedded in d e — 2 we find 
surprisingly, that although d — 2, de > 2 and d m in < 1> in contrast to regular lattices. 
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1. Introduction 

It has been realized in the last decades that a large number of complex systems are 
structured in the form of networks. The structures can be man-made like the World 
Wide Web and transportation or power grid networks or natural like protein and neural 
networks [?, 1-16, 19]. When studying the properties of these networks it is usually 
assumed that spatial constraints can be neglected. This assumption is certainly correct 
for networks like the World Wide Web (WWW) or the citation network where the real 
(Euclidean) distance does not play any role, but it may not be justified in networks 
where the Euclidean distance matters [20]. Typical examples of such networks include 
the Internet [6,11], airline networks [21,22], wireless communication networks [23], and 
social networks (like friendship and author networks) [24,25], which are all embedded in 
two-dimensional space (surface of the earth), as well as protein and neural networks [26], 
which are embedded in three dimensions. 

To model these networks, two network classes are of particular interest: Erdos-Renyi 
(ER) graphs [27,28] and Barabasi- Albert (BA) scale free networks [29]. In ER- networks, 
the distribution of the number k of links per node (degree-distribution) is Poissonian 
with a pronounced maximum at a certain /c-value, such that nearly each node is linked 
to the same number of nodes. In BA networks, the distribution follows a power law 
P(k) ~ k~ a , with a typically between 2 and 3. Here we focus on ER-type networks 
embedded in one- and two-dimensional space. We actually use a degree distribution that 
is close to a delta function (as the case in simple lattices). We found that the results are 
the same for both kinds of distributions. We follow Refs. [30-32] and assume that nodes 
are connected to each other with a probability p(r) ~ r~ 6 , where r is the Euclidean 
distance between the nodes. The choice of a power law for the distance distribution is 
supported from findings in the Internet, airline networks, human travel networks and 
other social networks [22,25,33]. Our model of embedding links of length r, chosen 
from Eq. (1), in a d e - dimensional lattice can be regarded as a generalization of the 
known Watts Strogatz (WS) model [1,3]. In the WS model links of any possible lengths 
with the same probability are added in the lattice system which corresponds to the case 
5 = of Eq. (1). Other methods for embedding networks in Euclidean space have been 
proposed in [34-37]. 

It has recently been shown that spatial constraints are important and may alter the 
dimension and therefore the topological properties of the networks (likethe dependence 
of the mean topological distance on the system size) as well as their robustness [30,31]. 
Here we are interested in studying how in these model networks the spatial constraints 
quantified by the distance exponent 5 modify the scaling relations between mass (number 
of nodes), Euclidean distance r and topological distance £. Our earlier study on ER 
networks embedded in a square lattice (with dimension d e = 2), indicate that by 
varying the exponent 5 one can actually change continuously the dimension d of the 
network, from d = oo for 5 < 2 to d = 2 for 5 > 4 [32]. In the present manuscript 
we present further extensive numerical simulations for d e = 2 that support this claim 
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as well as simulations in linear chains (d e = 1) that suggest analogous conclusions. In 
d e — 1 we find that for 5 < 1 the system behaves like an infinite dimensional network 
(as the original ER-network). When continuously increasing 5 the dimension becomes 
finite for 5 > 1 and approaches d = 1 for 5 > 2. Since the dimension of a system 
plays a critical role in many physical phenomena like diffusion, percolation and phase 
transition phenomena, our results are important for understanding and characterizing 
the properties of real world networks. 

Our manuscript is organized as follows. In Section 2, we discuss the characteristic 
distances in the spatially constrained networks. In Section 3 we describe the method 
to generate the spatial network models. In Section 4 we present our numerical results 
for the dimension d, for networks embedded in linear chains and in square lattices, that 
we obtain from the scaling relation of the mass M and the distance r. In Section 5 
we present our numerical results for the dimension d, that we obtain from the scaling 
relation of the probability of return to the origin P of a diffusing particle and its distance 
r. In Section 6, we discuss the scaling of the mass M and the Euclidean distance r with 
the topological distance t. The conclusions in Section 7 summarize our main results. 



5 =0.5d 8 = 1.5 d 8 = 2.5 d 

e e < 




Figure 1. Illustration of ER networks embedded in linear chains (top) and square 
lattices (bottom), for various distance exponents S. 



2. Characteristic distances 

First we estimate how the characteristic distances, in a network of L de nodes, depend on 
its linear size L, on 6 and on the embedding dimension d e . We normalize the distance 
distribution p(r) such that f. drr de_1 p(r) = 1, which yields 

(4 - 5)L-^-V r~ 5 ,5< d e 

P(r) = { (1) 
(5-d e )r~ d ,5>d e . 
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From p(r) we obtain r n = dr r de 1 r n p(r) and the related length scales r n = (r 11 ) 1 /". 
The maximum distance r max is determined by L de dr r de ~ l p(r) ~ 1. The results 
for r n and r max are 




L n /ln(L) ,5 = d e 




n ln(L) , <5 = d e + n 



de—5 



,5>d e + n 



v 



(i e +rt— 5 



and 
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,5 <2d e 



r 



max — 



L d e /(s-d e ) ^>2d e . 



(3) 



Accordingly, for 5 < d e all length scales (f n and r max ) are proporional to L, the spatial 
constraints are weak and the system can be regarded as an infinite dimensional system. 
On the other hand, for 5 > 2d e , f n /L and r max /L tend to zero in the asymptotic limit. 
In this case, we expect that the physical properties of the network are close to those 
of regular lattices of dimension d e . However, large finite size effects are expected for 5 
close to 2d e where r max /L decays only very slowly to zero. In the intermediate 5-regime 
d e < 5 < 2d e , r max scales as L, while f n /L tends to zero in the asymptotic limit. In this 
regime our simulation results (Chap. 4) suggest intermediate behavior represented by a 
dimension between d e and infinity that changes with 5. 

3. Generation of the networks 

The nodes of the network are located at the sites of a <i e -dimensional regular lattice, in 
our case a linear chain of length L (d e = 1) or a square lattice of size L x L (d e — 2). 
We assign to each node a fixed number k of links (in most cases, k — 4). Actually this 
network is a random regular (RR) network since all nodes have the same degree. It 
is expected (and we have also verified it numerically) that both networks, ER and RR 
with the same spatial constraints, are in the same universality class. 

To generate the spatially embedded networks, we use the following iterative 
algorithm: (i) We pick a node % randomly and choose, for one of its available ki links, a 
distance r (1 < r < L) from the given probability distribution p(r), Eq. (1). It is easy 
to see that the distance r can be obtained from random numbers < u < 1 chosen from 
the uniform distribution, by 



(ii) We consider all iV r nodes between distance r — Ar and r from node i, that are not 
yet connected to node i. Without loss of generality, we choose Ar = 1 for the linear 
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Figure 2. The distance distribution p(r) r de_1 for ER-networks embedded (a) in linear 
chains where d e — 1 and (b) in square lattices where d e — 2, when 5 — 0.5d e (circle), 
1.5d e (diamond), 2d e (triangle up), 2. hd e (triangle left) and k = 4. The numbers denote 
the slopes of P(r)r d "~ 1 , which are identical to the anticipated ones. For the same set 
of parameter as in (a) and (b), the panels (c) and (d) show the degree distribution 
p(k), which is = 1 for k = 4 and = otherwise. 



chain and Ar = 0.4 for the square lattice, (hi) We pick randomly one of these nodes j. 
If node j has at least one available link, we connect it with node i. If not, we do not 
connect it. Then we return to (i) and proceed with another randomly chosen node. At 
each step of the process, either 2 or zero links are added. For generating the network, 
we have typically performed 10 3 • L de trials. Due to the generation process, the nodes 
of the final network do not all have exactly the same degree, but the degree follows a 
narrow distribution with a mean k slightly below k — 4. Figure 1 illustrates the ER 
networks embedded in d e =l and d e =2 for 5 = 0.5d e , 1.5d e and 2.5d e . Figure 2 shows 
the actual narrow degree distribution as well as p(r) obtained in the simulations. 

4. The dimension of the networks 

For determining the dimensions of the spatially embedded networks, we follow the 
method developed by Daqing et al [32]. We use the fact that the mass M (number 
of nodes) of an object within an hypersphere of radius r scales with r as 

M ~ r d (5) 

where the exponent d represents the dimension of the network. When using this relation 
without taking into account the way the nodes are linked, one trivially and erroneously 
finds that the dimension of the network is identical to the dimension d e of the embedding 
space. 
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Figure 3. Illustration of the shells 5(0) (square), 5(1) (circle), 5(2) (triangle up), 
5(3) (diamond), 5(4) (triangle left), 5(5) (triangle down) and 5(6) (triangle right) for 
ER networks embedded in a square lattice with k = 4 (left panel) , and the mass M as 
function of I within this shells (right panel) . 

To properly take into account the connectivity, when considering the dimension of 
the network, we proceed as follows (see Fig 3): We choose a node as origin and determine 




r / r 



Figure 4. (a) The mass M as function of the relative distance r/r for ER networks 
embedded in linear chains with k = 4, for the system sizes TV = 10 5 , 10 6 andl0 7 with 
5 = 1.25, 1.5, 1.75, 2, 2.25, 2.5 (from left to right). The straight lines are best fits to 
the data that yield the dimension d of the network, (b) The same as panel (a) , but for 
N = 10 7 and k = 3,4and6. 

its nearest neighbors (referred to as shell 1) and their number S(l), the number of second 
nearest neighbors 5(2), and so on. Next we measure the mean Euclidean distance r{£) of 
the nodes in shell £ from the origin and determine the number of nodes M(£) = Yfi=i Sty) 
within shell £. To improve the statistics, we repeat the calculations for many origin nodes 
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and then average r(£) and M(£). To reduce finite size effects, we do not choose the origin 
nodes randomly in the underlying lattice, but from a region with radius L/10 around 
the central node. From the scaling relation between the average M and the average r, 
Eq. (5), we determine the dimension d of the network. 

Figure 4 shows the results for networks embedded in linear chains, for distance 
exponents 5 between 1.25 d e and 2.5 d e . In (a), we consider networks with k = 4 fixed 
and different system sizes (N = 10 5 , 10 6 and 10 7 ), while in (b) we consider networks 
with a fixed size N = 10 7 and various k values (k = 3,4,6). In both panels, we have 
plotted M as a function of r/f, where f = fi is the mean distance, see Eq. (2). 

Figure 4a shows that for 5 in the interesting regime between d e and 2d e , the curves 
for different N collapse nicely (For transparency, the curves (except 5 = 1.25) have been 
shifted along the x-axis by a factor of 10, 10 2 , 10 3 , 10 4 and 10 5 ). From the slopes of 
the straight lines, we obtain the dimension d = 4.64 (5 = 1.25), d = 2.12 (5 = 1.5) and 
d = 1.48 (5 = 1.75). For 5 > 2, the data starts to overshoot above some crossover value 
that increases with the system size and thus can be regarded as a finite size effect. To 
understand the reason for this crossover note that a node close to the boundary has 
a considerably higher probability to be linked with nodes closer to the center of the 
underlying lattice. As a consequence, for large shell numbers £, the mean Euclidean 
distance of the nodes from the origin node will be underestimated and thus the mass 
within large Euclidean distances overestimated. This effect is most pronounced in the 
linear chain, for intermediate 5-values, and gives rise to the overshooting of M(r) for S 
between 2 and 2.5, where d ~ d e . For 5 = 2.5 and iV = 10 7 , the total number of nodes 
in the spatially constrained network is well below N, since the network is separated into 
smaller clusters. For larger /c-values, this effect is less likely to appear. Figure 4b shows 
that the dimension of the networks does not depend on their average degree. The M(r) 
curves collapse for different k, and thus give rise to the same dimensions. This indicates 
the universality feature of the dimension. 

Figure 5 shows the corresponding results for networks embedded in square lattices 
(d e = 2), again for 6 exponents S between 1.25 d e and 2.5 d e , three network sizes 
(N = 9 • 10 4 , 10 6 and9 • 10 6 ), and three k values (k = 3,4,6). From the slopes of the 
straight lines we obtain d = 5.82 (5 = 2.5), d = 2.91 (5 = 3), and d = 2.21 (5 = 3.5). 
For S above 4, d is close to d e , as expected. The figure confirms that the finite size effects 
in d e = 2 are considerably less pronounced than in d e — 1, contrary to the intuition, 
since the linear size of the underlying embedding lattice is considerably higher in d e — 1 
than in d e = 2. As in d e — 1, the dimensions are independent of the mean degree of the 
networks. 

Figure 6 summarizes our results for the dimensions of the spatially embedded 
networks in the intermediate 5 regime between d e and 2d e , where the dimension is 
supposed to bridge the gap between d = oo for the unconstrained case 5 below d e and 
d = d e for the highly constrained case S above 2d e . The figure shows d — d e as a function 
of the relative distance exponent 8' = 5/d e for both considered lattices. The figure 
shows that in both cases, the curves approximately collapse to a single line which can 
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Figure 6. The difference between network dimension d and embedding dimension d e 
as a function of 8/d e for d e = 1 (circles) and d e = 2 (triangles). 



be represented by 

d - d * = c ^Ty 1<5 '< 2 (6) 

where c = 1.60. According to Eq. (6), d — d e diverges for 5' approaching the critical 
relative distance exponent 5' — 1. 
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5. The probability of return to the origin 

The network dimension plays an important role also in physical processes such as 
diffusion [38-40]. The probability Po(t) that a diffusing particle, after having traveled t 
steps, has returned to the origin, is related to the root mean square displacement r(t) 
of the particle by [32,40,41] 

P (t)~r(ty d . (7) 




r / r 



Figure 7. (a) The probability Pq that a diffusing particle is at its starting site, after 
travelling an average distance r, as a function of the relative distance r/f for ER 
networks embedded in linear chains with k = 4, for the system size N = 10 7 with 
5 = 1.25, 1.5, 1.75, 2, 2.25, 2.5 (from left to right). The straight lines are best fits to 
the data that yield the dimension d of the network, (b) The same as panel (a), but 
for ER networks embedded in a square lattice, the system size is N = 9 • 10 6 with 
S = 2.5, 3, 3.5, 4, 4.5, 5 (from left to right). Note that the values of d obtained here 
are almost the same as those obtained by direct measurements in Figs. 4 and 5. 

To derive Eq. (7) one assumes that the probability of the particle to be in any site in 
the volume V(t) = [r(t)] d is the same. As a consequence, Po(t) ~ 1/V(t), which leads to 
Eq. (7). Figure 7 shows Pq as a function of r/f in d e = 1 and 2, for the same 5- values as 
in Figs. 4 and 5. For convenience, we show only the results for the largest system size, 
N = 10 7 for d e = 1 and N = 9 ■ 10 6 for d e = 2. To obtain Po{t), we averaged, for each 
value of 5, over 10 4 diffusing particles and 50 network realizations. From the straight 
lines in the double-logarithmic presentations of Figure 7 we obtain the dimension of the 
networks, which are listed in the figure. The dimensions obtained in Figure 7 agree very 
well with those obtained by direct measurements in Figs. 4 and 5. 
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6. The topological dimension and the dimension of the shortest path 

In order to find how M scales with the Euclidean distance r, we determined in Sect. 4 
how M and r scale with the topological length £, and obtained the dimension d from 
M{£) ~ r(£) d . In this section, we discuss explicitely how M and r depend on I. 

It is well known that for regular lattices as well as for fractal structures, M and r 
scale with £ as power laws, 

M{£) ~ £ d ? (8 a) 

r{£) ~£ l l d ™^ (8b) 

where de is the topological ("chemical") dimension and d m ; n is the dimension of the 
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Figure 8. The mass M (left column) and the relative distance r/f (right column) 
as function of l a (i is the topological distance and a = S'(2 — 6')) for ER networks 
embedded in linear chains with k = 4, for the system sizes N = 10 5 , 10 6 andl0 7 with 
5 = 0.5, 1.25 and 1.75. The straight lines are best fits to the data with slopes A and B 
respectively. 



shortest path, see e.g., [42,43]. For regular lattices of dimension d e , di = d e and d min = 1. 
Thus we expect that for 5 > 2d e , the power law relations (8) hold. 

For 5 = the network has no spatial constraints and it is known that the mean 
topological distance (£) between 2 nodes on the network scales with the network size N 
as (£) ~ logiV [5]. This represents the small world nature of random graphs. Since N 
plays the role of the mass M of the network, it follows that M increases exponentially 
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with £, i.e. M(£) ~ exp(A£). We expect that this relation holds for 5 < d e where r max 
and r/f are both proportional to the linear scale L of the network, see Eqs. (2) and 
(3). Since for 5 > 2d e we expect power law relations (8), we conjecture that in the 
intermediate regime d e < 5 < 2d e , M(£) will increase slower than exponential and faster 
than a power law, via a stretched exponential, 

M{£) ~ exp(AT), d e <5< 2d e . (9) 

This function can bridge between the exponential behavior for 5 < d e and the power 
law for 5 > 2d e . For 5 approaching d e from above, a should approach 1, while for 5 
approaching 2d e from below, a should approach 0, consistent with a power law. The 
conjecture, Eq. (9) is supported by earlier numerical simulations [30] where it was 
found that in the intermediate regime, £ scales as (logiV)^, leading to a = 1/(3. On 
the basis of numerical simulations it was estimated [30], that a ~ 5(2 — 5) in d e — 1 
and a ~ 5(4 — 5) /A in d e = 2, which actually can be combined into a single equation, 
a = 5' (2 — 5'), when the relative distance exponent 5' = 5/d e is introduced. Thus our 
conjecture (9) becomes 

z M , 8' < 1 

M(£) ~ I , , (10) 
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Figure 10. The mass M (left column) and the relative distance r/r (right column) 
as a function of the topological distance £ for ER networks embedded in linear chains 
with k = 4, for the system sizes N = 10 5 , 10 6 andl0 7 with S = 2.0, 2.25and2.5. The 
straight lines are best fits to the data that yield the topological dimension di and the 
dimension of the shortest path d m i n ■ Note that the slopes below the crossover in (b) , 
(c), (e) and (f) of M and r vs I are the same. This yields d = 1 for all range of r as 
indeed seen in Fig. 4. 



where the prefactor A may depend on 5' and d e . To test this hypothesis, we have 
plotted, in Figs. 8, a, b, c (d e = 1) and Figs. 9 a, b, c (d e = 2), M(£) versus £ s '( 2 - s '\ 
in a semi-logarithmic fashion. The relative distance exponents 5' are 0.5, 1.25 and 1.75 
in both cases. The lattice sizes are the same as in Figs. 4 and 5. For 5' = 0.5 where 
the spatial constraints are irrelevant, we find logM ~ £, in agreement with (10). In the 
intermediate 5 regime 1 < 5' < 2 we find that logM ~ £ a , with a = 0.93 (5 = 1.25) and 
0.43 (5 = 1.75), also in agreement with (10). Accordingly, in the intermediate 5-regime, 
M(£) scales with the topological distance £ as a stretched exponential which serves as 
a "bridge" between the exponential behavior for 5 < d e and the anticipated power law 
behavior for S well above 2d e . 

Now the question arises how the power law in Eq. (5) that describes the scaling 
of M with r and the stretched exponential in Eq. (10) that describes the scaling of M 
with £, can be simultaneously satisfied. The only way to fulfill both equations is, that 
also r(£) is a stretched exponential with the same a in the intermediate regime i.e., 

r(£)~e mS ' {2 - S '\l<5' < 2, (11) 
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Figure 11. The same as Fig 10, but for ER networks embedded in a square lattice. 
The system sizes are N = 9 • 10 4 , 10 6 and 9 • 10 6 with 6 = 4.0, 4.5 and 5.0. The lines in 
(c) and (f) demonstrate for comparison slopes 2 and 1 respectively. 



and the ratio between the prefactors A and B should yield the dimension of the network. 
This is since M(£) ~ e M& (2 a ' = (e B£a (2 6 )s ) A / B ~ r d . Figs. 8 e, f and 9 e, f support the 
assumption (11). The prefactor B is obtained from the slopes of the straight lines in 
the figures and indeed the values of A/B are found to be identical to the values of the 
dimensions we obtained in the previous section. For 5 below d e (see Figs. 8d and 9d), 
r is independent of £ and M ~ e M (see Figs. 8a and 9a). 

For 5 > 2d e , we expect that M(£) and r(£) follow power laws, such that we can 
determine, from a double logarithmic plot, the chemical dimension di and the dimension 
of the shortest path, d m i n . Figures 10 and 11 show that this is the case. But surprizingly, 
for 5 > 2d e (but close to 2d e ), the values of d m i n and di do not agree with the values 
for the corresponding regular lattices. For 5 = 2d e , we obtain d% ~ 3.02 in d e — 1 
and di ~ 3.67 in d e = 2, significantly higher than the corresponding values di = 1 
and di = 2 in regular lattices. Furthermore, the dimension of the shortest path d min is 
considerably smaller than in regular lattices (d m i n = 1), d m i n = 1/2.65 = 0.38 in d e = 1 
and d min = 1/1-80 = 0.56 in d e = 2. Since M ~ £ dl ~ r dminde , the dimension d of the 
network for 5 > 2d e is simply d = d m i n d^ which yields d ~ 1.14 in d e — 1 and d ~ 2.04 
in d e = 2, in agreement with our results of Figs. 4-7. For 5 above 2d e we expect that 
di and d min accept the values of the corresponding regular lattices. Figure 10 shows 
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that this is indeed the case in d e = 1, with a pronounced crossover behavior for 5 = 2.25 
and 2.5. The crossover point decreases with increasing S. In d e = 2, in contrast, for 
5 = 2.25d e and 2.5d e the dimensions do not seem to reach their anticipated values d e = 2 
and d m i n — 1, even though d = 2 was obtained for both 5 values. Figure 11 does not 
suggest that this is a finite size effect since a bending down for larger system sizes cannot 
be seen similar to that in d e — 1. However, we cannot exclude the possibility that at 
very large system sizes that right now cannot be analyzed with the current state-of-the- 
art computers, there will be a crossover towards the anticipated values of de — 2 and 

dmin 1 • 

7. Summary 

In summary, we studied the effect of spatial constraints on complex networks where the 
length r of each link was taken from a power law distribution, Eq. (1), characterized by 
the exponent 8. Spatial constraints are relevant in all networks where distance matters, 
such as the Internet, power grid networks, and transportation networks, as well as in 
cellular phone networks and collaboration networks [6,11,20,23-25]. Our results suggest 
that for 5 below the embedding dimension d e , the dimension of the network is infinite 
as in the case of netwoks that are not embedded in space (represented by 5 = 0). For 
S between d e and 2d e , the dimension decreases monotonically, from d — oo to d — d e . 
Above 2d e , d = d e . We also studied how the mass M and the Euclidean distance r scale 
with the topological distance £. For 5 below d e , M increases exponentially with £, while 
r does not depend on I. For 5 between d e and 2d ei both the mass M and the Euclidean 
distance r increase with £ as a stretched exponential, with the same exponent a but 
different prefactors in the exponential. The ratio between these two prefactors yields the 
dimension of the embedded network. Exactly at 5 = 2d e , the exponent a becomes zero 
and M and r scale with £ as power laws, defining the exponents de and ci m i n , respectively 
similar to fractal structures [42,43]. While the dimension d is equal to d e , surprisingly 
de and d m i n do not have the values de = d e and d m i n = 1 that are expected for regular 
lattices. This effect seems to hold in d e = 2 also for 5 values somewhat greater than 2d e . 

Our results have been obtained for a nearly ^-functional degree distribution, but we 
argue that they are valid for any narrow degree distribution, like Possonian, Gaussian 
or exponential degree distribution since all those networks are expected to be in the 
same universality class. For power law degree distributions (scale free networks [29]), 
there may be differences for small values of 5, since it is known that nonembedded 
random graphs and scale free networks are in different universality classes [44,45]. In 
the relevant intermediate 5 regime (d e < 5 < 2d e ), we cannot exclude the possibility 
that the dimensions do not depend on the degree distribution. Indications are from 
measurements of the dimension of the airline network and the Internet [32]. Both are 
scale free networks, with 5 close to 3 (airline network) and 5 close to 2.6 (Internet). 
For the airline network, d is close to 3, while for the Internet, d is close to 4.5. These 
values are consistent with those obtained here for the ER-networks, with the same 5- 
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values. We have assumed a power law distribution, Eq. (1), for the link length. Other 
distributions are possible, for example an exponential distribution which holds for the 
power grid and ground transportation networks [20]. This case is equivalent to 5 = oo, 
since we have a finite length scale and thus the dimension d of the network is expected 
to be the same as the dimension of the embedding space d e . 

A power law distribution of Euclidean distances appears also in other physical 
systems where the present results may be relevant. For example, model systems where 
the interactions between particles decay as r~ s have been studied extensively for many 
years, for recent reviews on the statistical physics and dynamical properties of these 
systems, see [46,47]. Magnetic models on lattices with long range bonds whose lengths 
follow a power law distribution have also been studied, see e.g., [48]. In Levy flights 
and walks, the jump lengths follow a power law distribution. For reviews see [39,53,54]. 
Finally, it has been found that a power law distribution of link lengths with 5 = d e or 
d e + 1 (depending on the type of transport) is optimal for navigation [19,49-52]. 
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